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Sediment  Dispersion  Near 
Dredge  Pipeline  Discharge 
in  Laguna  Madre,  Texas 


PURPOSE:  Results  of  field  measurements  taken  near  a  pipeline  discharge  of  fine-grained  dredged 
material  into  shallow  water,  a  revised  numerical  model  formulation,  and  comparisons  between  field 
and  model  results  are  presented  herein.  A  numerical  model  for  the  simulation  of  underflow 
spreading  resulting  from  a  pipeline  disposal  was  previously  presented  in  TN-DOER-Nll  (Teeter 
2001).  That  model  is  updated  here,  guided  by  the  analysis  of  field  measurements.  An  analysis  of 
entrainment  of  underflow  material  into  the  water  column  under  the  action  of  wind- waves  is  also 
presented. 

BACKGROUND:  Initial  dispersion  of  dredged  material  after  pipeline  discharge  is  important  to 
deposit  area  and  susceptibility  to  erosion  or  resuspension.  Dispersion  consists  of  a  fluid-mud  gravity 
underflow  and,  possibly,  an  overlying  water-colunm  plume.  Fluid-mud  thickness,  concentration 
structure,  and  overlying  water-column  suspension  concentration  were  measured  in  shallow,  wind- 
exposed,  micro-tidal  Laguna  Madre,  TX,  within  about  500  m  of  where  a  dredge  pipeline  was 
discharging.  Depths  were  0.5  to  2  m  and  currents  were  weak.  The  dredged  material  had  a  median 
particle  size  of  4  to  5  pm.  Median  fluid-mud  thicknesses  were  0.45  m  of  which  the  top  60  percent 
was  interpreted  as  underflow  and  the  remainder  as  deposit.  Fluid-mud  concentration  at  the  upper 
surface  of  the  underflow  layer  was  about  3  dry-kg/m^  and  increased  exponentially  with  depth  to 
about  48  dry-kg/m^.  The  deposit  was  48  to  1 10  dry-kg/m^  solids. 

A  numerical  model  that  would  simulate  underflow  fluid-mud  spreading  resulting  from  a  pipeline 
discharge  was  developed  as  an  aid  in  the  diagnosis  and  interpretation  of  field  measurements.  The 
model  was  based  on  one-dimensional  equations  for  momentum  and  mass  conservation.  Model 
features  that  limited  entrainment  and  concentration  change  caused  by  deposition  were  incorporated 
as  indicated  by  field  observations.  A  plume  of  suspended  sediment  200  to  500  mg/L  above  ambient 
concentration  occurred  over  the  underflow  footprint,  with  resuspension  driven  by  wind- waves.  The 
development  of  a  point  model  of  the  water  column  overlying  a  fluid  mud  layer  was  based  on  a 
balance  between  entrainment  and  settling.  Settling  was  prescribed  on  the  basis  of  a  laboratory-de¬ 
veloped  functional  dependence  on  concentration.  Data  were  used  in  the  model  to  estimate 
coefficients  for  this  entrainment  process. 

INTRODUCTION:  Maintenance  dredging  of  the  Gulf  Intracoastal  Waterway  (GIWW)  annually 
involves  more  sediment  than  the  total  of  the  natural  sediment  inputs  to  Laguna  Madre  (Morton, 
Nava,  and  Arhelgar  2001).  Hydraulic  dredging  followed  by  open-water  disposal  through  a  pipeline 
is  the  most  common  dredging  method  in  shallow,  vegetated  Laguna  Madre,  TX.  There  are  concerns 
that  redispersed  dredged  material  in  Laguna  Madre  is  contributing  to  turbidity  and  is  limiting  light 
penetration  to  sea  grasses  over  the  long  term  (Onuf  1994).  The  extent  of  initial  dredged-material 
spreading  is  important  information  for  assessing  total  resuspension  and  predicting  possible  sediment 
impacts  outside  designated  disposal  areas. 
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Measurements  of  fluid-mud  thickness,  concentration  structure,  and  overlying  suspension  concen¬ 
tration  were  made  north  of  Port  Mansfield,  TX,  in  Laguna  Madre,  within  about  500  m  of  an  ongoing 
dredged-material  pipeline  discharge.  Such  information  is  scarce,  but  is  needed  to  improve  under¬ 
standing  of  the  behavior  of  such  discharges  and  subsequent  sediment  dispersion.  Field  measure¬ 
ments  were  compared  here  to  results  taken  from  a  near  field,  numerical  underflow-spreading  model 
for  the  purpose  of  improving  model  formulation.  Information  was  also  used  to  characterize  pipeline 
discharges  in  a  large-scale  numerical  sediment-transport  model  of  this  system  (Teeter  et  al.,  in 
preparation).  A  simple  water-column  point  model  was  also  used  to  estimate  coefficients  for  an 
entrainment  relationship  that  would  describe  the  flux  of  underflow  sediment  into  the  water  column. 
As  will  be  shown,  the  extent  of  the  underflow  affects  the  extent  of  any  surface  plume  of  suspended 
sediment  which  might  form  during  or  shortly  after  discharge  as  a  result  of  entraimuent  of  the 
underflow  into  the  overlying  water  column. 

The  183-km  section  of  Gulf  Intracoastal  Waterway  (GIW  W)  in  shallow  Laguna  Madre,  TX,  is 
located  between  Port  Isabel  and  Corpus  Christi  Bay,  or  geographic  coordinates  26.1°  to  27°  N,  and 
approximately  97.4°  W.  This  section  requires  about  1 .6  x  1 0^  m^  of  maintenance  dredging  annually. 
About  75  percent  of  the  1,500-km^  Laguna  Madre  is  vegetated,  and  seagrasses  are  sensitive  to 
underwater  light  conditions  (Dunton  1994).  A  series  of  studies  has  recently  been  undertaken  to 
evaluate  the  impact  of  dispersion  of  sediments  from  the  placement  areas  on  sea  grass  areas  in  Laguna 
Madre  (Brown  and  Kraus  1997;  Militello,  Kraus,  and  Kite  1997;  Burd  and  Dunton  2000;  Morton, 
Nava,  and  Arhelgar  2001 ;  and  Teeter  et  al.,  in  preparation). 

BACKGROUND  ON  PIPELINE  DISCHARGE  UNDERFLOWS:  In  shallow  water,  dredged 
sediment  particles  reach  the  bottom  soon  after  pipeline  discharge,  settling  within  a  short  distance 
fi-om  the  discharge  point.  Sediments  form  layers  of  fluid  mud  at  the  bed,  which  flow  away  from 
the  point  of  discharge,  the  extent  of  the  flow  depending  on  bottom  slope,  ambient  currents,  and  their 
initial  discharge  trajectory.  Some  general  characteristics  of  underflows  were  presented  in  TN- 
DOER-Nll  (Teeter  2001). 

The  approach  channel  to  the  Chesapeake  and  Deleware  Canal  in  Upper  Chesapeake  Bay  was 
hydraulically  dredged  in  1988.  About  5.2  x  10^  m^  of  clayey  silt  sediment  were  pumped  and 
deposited  in  areas  D,  E,  and  F  near  Pooles  Island.  The  movement  of  sediment  w'as  downslope  after 
discharge.  A  broad  continuous  layer  formed  about  3  km  long  and  1 .5  km  wide.  The  maximum 
deposit  thickness  was  1.5  m.  Sediment  consolidated  to  a  density  of  1,130  kg/m^  or  greater  within 
several  weeks.  Dewatering  and  compaction  accounted  for  5  percent  deposit-volume  reduction 
in  5  months.  Another  5  percent  reduction  occurred  during  the  discharge  period.  The  remaining 
22  percent  of  the  32  percent  total  reduction  was  from  redistribution  by  resuspension  and  transport 
(Panageotou  and  Halka  1 990). 

Near  Pooles  Island,  5.2  x  10^  m^  were  hydraulically  dredged  from  the  nearby  channel  in  1991  and 
placed  in  Areas  D  and  E.  Sites  were  4.5  to  8  m  deep.  Sediments  deposited  in  a  natural  trough  and 
constructed  trenches.  The  sediment  remained  in  the  deep,  trough  area.  The  volume  of  the  deposit 
was  1.04  X  10^  m^  with  maximum  thickness  of  3  m.  Sediment  were  clayey  silt  with  minor  sand. 
Bulking  factor  between  in-place  and  deposited  volumes  was  about  1 .75.  One  year  later  the  deposit 
was  4.4  X  10^  m^  (58  percent  reduction).  Four-fifths  of  the  reduction  was  attributed  to  dewatering, 
one-fifth  to  erosion  (Panageotou  and  Halka  1 994). 
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FIELD  OBSERVATIONS:  Field  experiments  were  carried  out  to  take  advantage  of  dredging 
conducted  in  February  2000  to  remove  a  1 .8-m  layer  of  material  deposited  in  the  GIWW  as  a  result 
of  a  hurricane  the  previous  year.  Dredge  J.N.  Fisher  discharged  into  open- water  disposal  sites 
through  a  50. 8-cm-diam  pipeline,  using  a  1 ,5 00-kW  (2,000-hp)  pump.  The  dredging  rate  was  about 
1,100  m^/hr,  and,  based  on  the  solids  content  of  the  channel  material,  the  sediment  discharge  rate 
was  about  50  dry-kg/sec. 

Fluid-mud  thicknesses,  or  heights,  and  densities  were  measured  on  2  days  while  pipeline  discharge 
was  occurring.  Locations  for  the  discharges  are  shown  in  Figure  1.  A  special  push-tube  sampler 
allowed  for  fluid-mud  density  determination  within  only  a  few  minutes  of  sampling.  Samples  were 
collected  for  analysis  of  the  fluid-mud  concentration  to  supplement  the  field-density  measurements. 
Fluid-mud  particle-size  distribution  and  ambient  water  column  suspended-sediment  concentrations 
were  also  measured.  A  composite  sample  was  used  in  the  laboratory  to  determine  velocities  in  the 
hindered  settling  range. 

Field  Methods.  A  5.8-m-long  flat-bottomed  boat  with  a  propeller  tunnel  to  minimize  draft  was 
used  for  sampling.  A  Starlink  ®  Differential  Global  Positioning  System  was  used  to  locate  stations 
to  within  ±2  m,  and  an  HP  PalmPC  ®  was  used  to  log  positions  in  the  field.  Water-column  samples 
were  collected  with  a  submersible  Rule  ®  electric  pump  and  1 .5-cm-diam  hose.  Water  samples  were 
collected  at  middepth  and  0.3-m  depth  and  stored  in  225-ml  plastic  bottles. 

Fluid  mud  was  sampled  with  a  push  corer  with  a  clear  3.6-cm-diam  core  tube  and  a  total  length  of 
about  3  m.  During  the  first  sampling  day,  it  was  found  that  the  in-line  check-valve  developed  too 
much  back  pressure,  resulting  in  significant  errors  in  underflow  sampling.  For  the  second  sampling 
day,  fluid-mud  samples  were  collected  with  a  low  back-pressure  push-core  sampler  specially 
fabricated  from  parts  of  a  WILDCO  ®  corer.  That  sampler  can  be  seen  in  Figure  2.  Only  the  fluid 
mud  measurements  from  the  second  day  are  reported. 

The  boat  was  brought  to  a  new  location,  and  the  anchor  was  set.  A  couple  of  minutes  were  allowed 
for  the  boat  to  swing  to  and  for  the  position  to  be  logged.  The  corer  was  pushed  vertically  downward 
by  hand  until  it  encoimtered  firm  bottom.  A  trip  line  was  then  pulled  to  seal  the  top  of  the  sampling 
tube.  The  corer  captured  ambient  water  column,  fluid  mud  (if  present),  and  a  short  plug  of  the 
underlying  bottom  material.  (The  bottom  material  contained  an  appreciable  sand  fraction  not 
present  in  the  dredged  material  and  had  a  bulk  density  of  roughly  1,500  kg/m^).  The  vertical 
alignment  of  the  core  tube  was  maintained  as  it  was  lifted  to  the  deck  and  a  piston  push-rod  was 
inserted  into  the  lower  end  of  the  core  tube  (below  the  sediment  plug).  After  the  core  tube  was 
unscrewed  from  the  remainder  of  the  sampler,  the  piston  rod  was  pushed  upward  to  expel  the  sample 
from  the  end  of  the  tube.  By  incrementally  extruding  the  sample  from  the  end  of  the  core  tube, 
scientists  could  take  measurements  and  subsamples  over  the  vertical  dimension  of  the  fluid  mud. 
Density  measurements  were  made  in  the  field  with  a  PARR  ®  DM  A3  5  vibrating-tube  densitometer 
(precision  of  1  kg/m^).  A  short  length  of  2-mm-diam  tubing  was  inserted  2.5  cm  into  the  end  of 
the  core  tube,  and  a  5  to  10  cm^  sample  was  drawn  through  the  densitometer.  Field  density 
measurements  were  made  in  duplicate  and  averaged. 

Laboratory  Methods.  Laboratory  bulk  wet  density  determinations  were  made  with  the  use  of 
25  -cm^  wide-mouth  pycnometers.  Pycnometers  were  weighed  after  being  mostly  filled  with  sample 
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Figure  1 .  Vicinity  sketch  of  Lower  Laguna  Madre  north  of  Port  Manfield,  TX,  with  depth  contours  and 
discharge  locations  (coordinates  are  State  Plane  NAD27,  Texas  South,  in  meters) 
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Figure  2.  Fluid  mud  sampler  on  deck  in  open  position 
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and  then  carefully  topped  with  distilled  water.  Bulk  wet  density  was  calculated  from  this  informa¬ 
tion  and  known  characteristics  of  the  pycnometers.  Pore-fluid  density  was  estimated  on  the  basis 
of  the  salinity  determined  on  suspended  samples,  allowing  the  calculation  of  sample  solids  content 
from  bulk  wet  density  (assuming  a  solids  density  of  2,650  kg/m^). 

Total  suspended  material  (TSM)  was  determined  by  a  gravimetric  method  for  nonfilterable  solids 
with  preweighted  Nuclepore  ®  0.45  pm  pore  diam,  polycarbonate  filters.  After  being  used  to  filter 
a  known  volume,  filters  were  rinsed  with  distilled  water,  and  dried  1  hr  at  90  °C  and  then  reweighed. 
Particle-size  distribution  was  measured  with  a  Coulter  LSIOOQ  ®  laser  scattering  instrument. 
Samples  were  first  oxidized  with  Clorox  ®  to  remove  organics  and  then  were  dispersed  with  sodium 
carbonate/bicarbonate.  Three  oxidation  steps  and  three  dispersion  steps  were  performed  before 
samples  were  processed  through  the  Coulter  instrument  to  determine  particle  size.  The  Coulter  has 
128  geometrically  spaced  channels,  or  bins,  for  sizing. 

Settling  velocities  in  the  hindered-settling  concentration  range  were  measured  on  left-over  sample 
that  had  been  composited  to  make  a  slurry.  The  slurry  had  a  bulk  density  of  1,109.5  kg/m^ 
pore-fluid  density  of  1 ,025.7  kg/m^  (37.3  ppt),  and  solids  content  of  1 36.7  dry-kg/ml  Sample  was 
incrementally  added  to  a  2-L  glass,  graduated  cylinder  which  was  7.74  cm  in  diam  and  42.5  cm 
high  at  the  2-L  level.  Six  tests  with  concentrations  of  6.8  to  66.3  dry-kg/m^  were  made  at  23  °C. 
After  the  sample  was  mixed  in  the  cylinder,  height  of  the  interface  between  the  suspension  and  the 
clear  layer  that  formed  was  observed  over  time.  The  duration  of  the  lowest  initial-concentration  test 
was  about  1  hr.  During  other  tests,  frequent  measurements  were  collected  over  100  to  240  min; 
these  tests  lasted  a  total  of  1 , 1 00  to  1 ,450  min.  Final  data  points  allowed  for  estimation  of  average 
density  after  about  1  day  of  settling  time. 

Linear  regressions  were  fit  to  the  data  for  the  period  when  the  interface  descended  linearly  (/?  =  3 
to  20,  E? = 0.944  to  0.999,  standard  error  on  slope  =  1 .82  to  0.025  mm/min)  to  determine  the  hindered 
settling  velocity  (W^)  at  initial  test  concentration  C.  Tests  with  the  lowest  two  concentrations  were 
repeated,  and  data  sets  were  combined  in  the  regression  analysis.  Finally,  and  initial  concentra¬ 
tion  from  the  six  tests  were  combined  and  fit  to  an  empirical  equation  for  hindered  settling 
dependence  on  concentration: 

=  Wh^  {1  -k  Cy  ,  C  >  hindered-settling  threshold  (1) 

where  the  hindered  settling  threshold  is  usually  in  the  range  of  1  to  10  dry-kg/m^. 

Field  and  Laboratory  Results.  All  settling  tests  were  in  the  hindered  settling  concentration 
range.  Settling  rates  decreased  about  two  orders  of  magnitude  over  the  concentration  range  tested. 
Data  greater  than  6.8  dry-kg/m^  fit  Equation  1  well  with  the  reference  hindered  settling  velocity 
Wh^  =  0.5  mm/sec,  coefficient  k  =  0.005  m^/kg,  and  the  exponent «  =  1 1 .  Settling  test  results  are 
plotted  in  Figure  3  along  with  results  from  the  low-concentration  settling  tests  performed  by  Teeter 
et  al.  (in  preparation),  using  Laguna  Madre  GIWW  sediments  collected  about  3  km  north  of  Port 
Isabel.  The  mean  depth-average  concentration  at  the  end  of  the  settling  tests  (about  20  hr)  was 
1 15.5  dry-kg/m^  (with  one  high  outlier  of  148  dry-kg/m^  removed,  n  =  4,  95  percent  confidence 
interval  1 12.0  to  1 18.9  dry-kg/m^). 
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Figure  3.  Hindered  settling  test  results  (right)  with  the  (dashed)  fit  to  the  data  described  in  the  text,  and 
low-concentration  settling  results  (left)  from  Teeter  et  al.  (in  preparation) 

February  10.  When  sampling  began  at  1000  Central  Standard  Time  (CST)  on  10  February,  south 
winds  were  9  m/sec,  making  sampling  conditions  difficult.  Currents  were  weak  and  toward  the 
south.  The  pipeline  discharge  was  located  at  coordinates  26°  37.4358’  N  and  97°  24.8643’  W,  about 
3.7  km  north  of  the  entrance  to  Port  Mansfield  at  the  disposal  area  designated  PA  218  (Figure  1).  A 
plot  of  the  station  locations,  CST,  and  depths  in  a  local  horizontal  coordinate  system  are  shown  in 
Figure  4.  The  pipeline  was  located  20  m  west  of  the  1016  CST  station  and  changed  only  slightly 
as  the  dredge  moved.  Water-column  samples  taken  at  0.3  m  depth  had  TSM  levels  (mean  =  211 
mg/L,  95  percent  confidence  interval  51  to  370  mg/L,  n  =  12)  equivalent  to  those  from  middepth 
(mean  =  199  mg/L,  95  percent  confidence  internal  55  to  344  mg/L,  n  =  1 1).  Both  sampling  depths 
showed  highly  variable  TSM.  Depth-averaged  TSM  are  shovm  in  Figure  5.  High-suspension 
concentrations  were  measured  both  north  and  south  of  the  pipeline  discharge.  Depth-average  TSM 
values  of  549,  557  and  572  mg/L  were  obtained  near  and  downdrift  within  about  660  m  of  the 
discharge.  Stations  taken  at  1 107, 1 126,  and  1255  CST  had  what  was  apparently  background  level 
TSM  ranging  from  about  1 00  to  1 20  mg/L.  One  station  ( 1 200  CST)  taken  upstream  of  the  discharge 
had  262  mg/L  TSM,  possibly  as  the  result  of  local  resuspension.  ^ 

A  photograph  taken  from  1,900-m  altitude  above  the  dredging  operation  at  1019  CST  is  shown  in 
Figure  6.  The  pipeline  length  was  about  450  m,  and  the  discharge  flowed  to  the  east  of  the  channel. 
Ambient  and  dredged-material  plumes  at  both  the  dredging  and  disposal  sites  are  shown  in  the 
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Figure  4.  Times/depths  (CST/m)  for  stations  taken  1 0  February  2000 
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Figure  5.  Depth-averaged  TSM  (mg/L)  collected  10  February  2000 
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photograph.  The  aerial  photograph  shows  the  area  directly  south  of  the  discharge  to  be  the  heaviest 
visible  plume.  A  plume  emanating  directly  from  the  discharge  had  a  blue  coloration,  while  other 
plume  areas  were  milky. 

February  16.  Winds  were  from  the  south.  Waves  were  0.30  cm  or  less.  Currents  were  weak  and 
moved  toward  the  north  (~2  cm/sec).  The  pipeline  discharge  was  located  at  coordinates  26° 
44.9752’  N  and  97°  27.3349’  W,  about  3.3  km  south  of  the  entrance  to  the  Land  Cut,  at  PA  213 
(Figure  1).  The  discharge  was  to  the  east  of  the  channel  onto  a  mound  with  about  a  0.3-m  water 
depth.  Dredging  records  indicated  that  the  previous  discharge  location  was  360  m  north  in  PA  213. 
Discharge  started  there  at  about  2200  CST  on  15  February;  the  discharge  that  was  sampled  began 
at  about  0800  CST  on  16  February.  During  the  sampling  period,  12,000  to  15,000  m^  of  dredged 
material  was  discharged  at  these  two  sites. 

Station  times  and  depths  are  shown  in  Figure  7.  A  turbulent  surface  flow  formed  in  the  vicinity  of 
the  discharge  jet  and  extended  into  deeper  water.  A  photograph  of  the  surface  jet  and  flow  is  shown 
in  Figure  8.  A  plunge  line  could  be  clearly  seen  in  the  field  at  a  water  depth  of  about  1  m,  and  an 
underflow  moved  toward  the  deeper  water  to  the  east-southeast.  Samples  taken  at  0915  CST  were 
within  the  turbulent  surface  flow,  and  two  field  measurements  and  two  pump  samples  indicated  that 
the  turbulent  surface  flow  averaged  17  dry-kg/m^. 
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Figure  7.  Times/depths  (CST/m)  for  stations  taken  1 6  February  2000 
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Plunge  Line 


Figure  8.  Pipeline  discharge  into  approximately  0.5-m  water  depth  and  resulting  turbulent  surface  flow 


Fluid  mud  formed  a  sharp  interface  with  the  ambient  suspension,  and  its  thickness  was  easily 
measured  through  the  clear  3.6-cm-diam  core  tube.  Fluid-mud  profiles  are  presented  in  Figure  9. 
Because  of  concerns  about  settling  effects  and  time  constraints,  few  samples  near  the  upper 
underflow  interface  were  made.  The  fluid-mud  layer  was  highly  stratified  in  the  vertical.  Gradients 
indicated  that  concentrations  at  the  upper  interface  were  low.  Several  measurements  indicated 
minimum  underflow  concentrations  of  3  to  5  dry-kg/m^.  However,  the  upper  surface  of  the 
underflow  was  a  distinct,  sharp  interface,  indicating  a  concentration  jump  associated  with  the 
maximum  flux  of  suspended  material  (Teeter  1986).  Concentrations  in  the  underflow  were  therefore 
above  the  concentration  at  which  the  maximum  settling  flux  occurred.  The  maximum  settling  flux 
apparently  occurred  between  the  settling  flux  at  6.8  dry-kg/m^  (0.0051  kg/m^/sec)  determined  in 
the  hindered  settling  tests  and  the  previous  1  dry-kg/m^  low-concentration  test  (0.0006  kg/m^/sec). 
Thus,  3  dry-kg/m^  was  estimated  to  be  the  minimum  or  underflow  interface  concentration.  Points 
were  added  to  profiles  as  shown  in  Figure  9  at  the  measured  interface  locations  and  the  assumed 
3  dry-kg/m^. 

Fluid-mud  layers  consisted  of  underflow  and  deposit,  as  interpreted  by  the  following  information. 
Near  the  top  of  the  fluid-mud  layers,  concentrations  increased  exponentially  with  depth  and  in 
approximately  straight  lines  when  plotted  on  semi-log  axes  as  in  Figure  9.  This  distribution  would 
be  expected  for  a  turbulent  flow  with  a  particle  Peclet  number  (Pe  =  W^h/K^  where  h  is  the  underflow 
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Figure  9.  Vertical  fluid  mud  profiles  taken  1 6  February  2000 


thickness  and  is  the  layer-average  turbulent  diffiisivity)  greater  than  about  one  (Teeter  1986). 
Another  steeper  gradient  was  evident  below  many  of  these  exponential  layers,  and,  taking  all 
measurements  together,  the  statistical  distribution  of  fluid-mud  concentrations  had  an  inflection  at 
about  50  dry-kg/m^.  Previous  laboratory  experiments  on  sediment  from  nearby  Corpus  Christi  Bay 
(76  percent  clay  and  21  percent  silt)  indicated  that  the  mean  concentration  of  newly-deposited 
material  was  46  dry-kg/m^  (Teeter  1986).  Therefore,  the  lower  portion  of  these  fluid-mud  layers 
was  interpreted  as  deposit  from  the  underflow,  while  the  upper  layer,  with  concentrations  of  about 
3  to  48  dry-kg/m^,  was  interpreted  as  the  underflow.  Normalized  underflow  concentration  profiles 
were  similar,  as  can  be  seen  in  Figure  10,  and  their  exponential  shapes  suggest  some  degree  of 
vertical  mixing,  consistent  with  a  flow  with  some  turbulence. 

All  measured  fluid-mud  thicknesses  and  a  rough  interpretation  of  the  underflow  footprint  extent  are 
shown  in  Figure  11.  The  0915  CST  samples  were  assigned  0.0  fluid-mud  thickness  in  this  figure 
because  they  were  located  within  the  turbulent  surface  flow.  It  appeared  that  the  underflow  footprint 
formed  by  the  discharge  ongoing  during  sampling  overlapped  that  formed  at  the  previous  discharge 
location.  Underflow  mean  concentrations  C,  thickness  h,  and  deposit  thickness  delbed  were 
calculated  on  the  basis  of  a  50-point  interpolation  over  the  fluid-mud  profiles  as  shown  in  Figure 
12  for  these  stations,  along  with  the  water  column  depth  above  the  underflow 

Surface  TSM  levels  (mean  =  258  mg/L,  95  percent  confidence  interval  1 14  to  402  mg/L,  n  =  17) 
and  middepth  levels  (mean  =  262  mg/L,  95  percent  confidence  interval  -205  to  728  mg/L,  n  =  8) 
were  equivalent  again  this  day.  A  plot  of  depth-mean  TSM  values  is  shown  in  Figure  13,  along 
with  an  interpreted  underflow  footprint  extent.  As  can  be  seen,  the  highest  suspended  concentrations 
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Figure  10.  Interpreted  and  interpolated  underflow  profiles  normalized  by  underflow  thicknesses  and  mean 


concentrations 


Figure  1 1 .  Fluid  mud  thicknesses  collected  1 6  February  2000 
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Underflow  Mean  C,  kg/cu  m  Underflow  Thickness  h,  m 


300  500  700  900  300  500  700  900 

EASTING,  m  EASTING,  m 

Deposit  Thickness,  m  Depth  Over  Deposit,  m 


300  500  700  900  300  500  700  900 

EASTING,  m  EASTING,  m 


Figure  12.  Interpolated  underflow  mean  concentrations,  thicknesses,  deposit  thickness,  and  water 
column  depth  over  underflow  for  mud  profile  stations,  16  February  2000 


Figure  13.  Depth-averaged  TSM  (mg/L)  collected  16  February  2000 
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occurred  over  the  underflow,  whether  upstream  or  downstream  from  the  discharge  point.  The 
implication  is  that  wind-waves  were  acting  to  entrain  material  from  the  active  underflow  and  from 
the  previous  underflow  into  the  water  column.  Since  water-column  advection  was  minimal,  the 
resulting  TSM  plume  did  not  appear  to  extend  much  beyond  the  footprint  extent  (and  vice  versa). 

Comparisons  between  field  measurements  and  laboratory  densities  indicated  that  many  field 
samples  had  sampling  errors  and  were  biased  toward  lower  density.  These  samples  took  much 
longer  to  obtain  than  the  field  density  measurements  did,  and  settling  may  have  caused  the  bias. 
Results  from  the  field  densitometer  and  pycnometers  agreed  well  on  samples  tested  in  the  laboratory . 
Therefore,  density  results  from  samples  were  not  included  in  previous  figures.  Table  1  summarizes 
laboratory  analyses  of  particle-size  characteristics  on  fluid-mud  point  samples.  In  that  table,  H  is 
water  depth  to  the  original  bed,  z,  is  the  distance  up  from  the  bottom,  and  D50  is  the  median  dispersed 
grain  diameter. 


Table  1 

Summary  of  Fluid-Mud  Sediment  Characteristics  from  Near  Pipeline  Discharge 

Time,  CST 

H,  m 

Zi,  m 

D50.  pm 

%  <4  pm 

%  <16  pm 

915 

0.9 

0.6 

4.4 

47 

84 

915 

0.9 

0-0.15 

5.1 

43 

79 

933 

1.8 

0-0.15 

5.1 

43 

78 

1007 

1.5 

0.15 

4.2 

49 

86 

1007 

1.5 

0-0.15 

4.1 

50 

85 

1057 

1.9 

0.3 

3.9 

50 

87 

1057 

1.9 

0-0.15 

4.3 

49 

84 

1110 

1.8 

0-0.15 

4.4 

47 

85 

1122 

1.1 

0-0.15 

4.2 

49 

86 

1136 

2 

0-0.15 

4.1 

49 

86 

1216 

2 

0-0.15 

4.1 

50 

86 

1233 

1.7 

0-0.15 

4 

50 

88 

1308 

2.1 

0-0.15 

4.2 

50 

86 

UNDERFLOW  SPREADING  PROCESSES:  To  better  imderstand  the  spreading  of  the  fluid- 
mud  underflow,  a  mathematical  description  of  underflow  processes  was  developed.  Unfortunately, 
there  are  no  analytic  solutions  for  the  case  of  a  particle-driven  gravity  flow  which  is  entraining  and 
depositing  material,  so  a  numerical  solution  was  developed.  Important  model  features  were  guided 
by  field  observations  with  special  attention  to  entrainment  and  settling,  either  of  which  can 
appreciably  reduce  underflow  concentration. 

Model  Description.  A  model  was  constructed  to  compute  total  flow  or  discharge  (Q),  sediment 
flux  (C0,  breadth  (5),  and  height  {h)  along  the  length  (x)  of  an  underflow  by  numerically  integrating 
a  set  of  governing  equations  downslope  in  the  direction  of  the  underflow.  The  development  of  both 
model  equations  and  assumptions  were  guided  by  the  field  observations.  The  behavior  of  mobile 
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fluid-mud  layers  is  an  active  area  of  research,  and  more  studies  will  be  needed  before  such 
computations  can  be  made  with  confidence. 

The  bed  was  assumed  to  be  planar  with  an  arbitrary  slope  which  was  allowed  to  vary  in  the 
downslope  direction.  The  underflow  was  considered  quasisteady  or  steady  over  a  short  duration  of 
time.  Thus,  time  derivatives  were  ignored  in  the  governing  equations.  However,  time  of  travel  to 
every  discrete  location  along  the  underflow  trajectory  was  calculated  by  integrating  the  underflow- 
velocity.  After  the  first  full  numerical  integration  sweep  from  the  transition  (beginning)  to  the  end 
of  the  underflow,  subsequent  sweeps  were  made  at  discrete  time  intervals  and  included  updated  bed 
elevations  based  on  the  cumulative  deposit  thickness  from  preceding  sweeps.  Sweeps  were  made 
at  intervals  of  1,800  sec  so  that  per-sweep  deposit-thickness  changes  were  small.  Thus,  the  model 
made  many  sweeps  over  the  underflow  domain  and  duration  of  the  discharge  to  update  deposit 
thickness  delbed  and  bed  slope.  Other  variables  were  calculated  from  the  basic  state  variables  and 
used  to  solve  auxiliary  equations  for  entrainment,  bed  friction  coefficient,  depositional  flux,  and 
lateral  spreading.  These  variables  include  underflow  concentration  (Q,  velocity  (60,  deposit 
thickness,  Reynolds  and  Richardson  numbers. 

Equations.  Governing  layer-averaged  equations  are  based  on  mass  and  momentum  conservation. 
Governing  equations  for  steady,  one-dimensional  gravity  flow  momentum  have  been  developed  by 
Ellison  and  Turner  (1959),  Parker,  Fukushima,  and  Pantin  (1 986),  and  Van  Kessel  and  Kranenburg 
(1996).  Ellison  and  Turner  (1959)  ignored  Ap/pj  terms  except  when  combined  with  g  where^  = 
the  acceleration  of  gravity,  Ap  =  the  density  difference  between  the  underflow  and  the  overlying 
ambient  suspension  (i.e.,  p  -  p,  where  P|  =  constant ),  p  =  the  layer  density  (p  =  P|  +  Ap).  This  is 
equivalent  to  ignoring  them  in  inertia  terms  of  the  equations  of  motion,  and  thus  the  fluid  is  treated 
as  having  uniform  mass,  variable  weight,  and  conservative  buoyancy.  Parker,  Fukushima,  and 
Pantin  (1986)  derived  equations  explicitly  for  particle-driven  gravity  currents  where  deposition  and 
erosion  can  occur  and  can  affect  layer  buoyancy.  They  began  with  the  three-dimensional  Navier- 
Stokes  equation.  However,  they  also  assume  uniform  mass  density.  Van  Kessel  and  Kranenbrng 
(1996)  cite  Ellision  and  Turner  (1959)  as  developers  of  an  appropriate  layer-averaged  integral 
momentum  equation  but  made  a  key  change  by  keeping  the  density  term  inside  the  momentum 
derivative.  In  this  way  mass  density  and  mass  flux  are  not  assumed  constant  in  the  x-direction.  The 
following  paragraphs  derive  a  form  of  their  equation  for  a  wide,  variable-width  flow  of  rectangular 
cross  section. 

The  momentum  equation  was  developed  considering  the  forces  acting  on  the  boundaries  of  an 
elemental  flow  volume.  Some  definitions  are  given  in  Figure  14.  The  momentum  of  a  flow  volume 
with  width  B  between  x  and  x  +  cfr  and  planes  normal  to  the  x-axis  and  parallel  to  the  z-axis  is  the 
mass  times  its  velocity.  Let  p  be  density  of  the  flow.  Then  p5«  dt  is  the  mass  passing  through  the 
upstream  plane  per  unit  depth,  where  u  is  the  mean  flow  across  the  plane.  Since  the  mass  moves  at 
u,  the  time-averaged  momentum  is  pBu^  per  unit  depth. 

Consider  the  forces  acting  on  the  boundaries  of  the  flow  volume.  They  include  the  pressure  forces 
(in  excess  of  the  hydrostatic  pressure)  acting  on  the  downstream  and  upstream  planes,  reactive 
shear-stresses  acting  along  the  bottom  and  the  top  (T|,  +  Xj  =  t  where  =  bed  shear  stress,  t,  —  shear 
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Figure  14.  Definitions  for  some  model  variables 


stress  at  the  top  of  the  underflow,  and  x  is  the  total)  ,  and  a  body-force  accelerating  the  volume 
downslope.  Then  the  momentum  balance  is 


A 

dx 


I  p  dx  =  - -  g  cos 0  J  A  Bxdz  +  g  sin  0  I  Ap  Bdz 


(2) 


where  the  left-hand  side  term  is  the  change  in  momentum  along  the  slope,  and  the  right-hand  side 
terms  are  (a)  the  total  shear  stress  resisting  the  flow;  (b)  the  derivative  of  the  pressure  along  the  flow 
resulting  from  changing  depth  and  density;  and  (c)  the  gravity  force  acting  to  accelerate  the  layer. 

Assume  h(x),  Ap(x),  B(x),  and  U(x)  vary  along  the  flow.  Define  variables  for  mass  and  momentum: 
pBUh  =  j  pBu  dz  and  pBU^h  =  j  pBu^  dz.  Shape  factors  S'!  2  are  defined  during  integration  of  pres¬ 
sure  and  gravity  terms  taking  into  account  z-direction  variations  in  density  p(x)  such  that 


SxApBh^l2  =  ^ApBzdz 
$2  Ap  B/i  =  J  ApBdz 


(3a,b) 


The  excess  density  is  integrated  and  the  hydrostatic  pressure  does  not  appear.  Thus,  vertical 
integration  is  done  under  the  assumption  that  shape  factors  do  not  vary  in  the  x-direction.  The 
momentum  equation  is  then  cast  as 


d(^BU'^}^ 


dx 


=  -Bx 


5igcos0 


l^pBh^^ 


dx 


+  S2gApBhsmQ 


(4) 


where  0  =  the  angle  of  the  bed  from  the  horizontal. 
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Van  Kessel  and  Kranenburg  ( 1 996)  assumed  the  shape  factors  were  unity  for  laminar  flow  and  used 
a  form  of  Equation  4  with  5  =  1  to  analyze  turbulent  flows.  Parker,  Fukushima,  and  Pantin  (1986) 
evaluated  shape  factors  using  experimental  data  and  found  the  top  hat  assumption  for  all  vertical 
profiles  (including  momentum)  to  be  a  good  approximation.  Ellison  and  Turner  (1959)  found 
-Si  =  0.2  to  0.3  and  S2  ~  0.6  to  0.9  in  their  experiments.  Density  profiles  from  Laguna  Madre  were 
integrated  after  exponential  interpolation  at  1-cm  increments,  between  measurement  points  («  =  19 
to  47).  Profiles  are  shown  in  Figure  9  and  dimensionless  profiles  are  shown  in  Figure  10.  The 
trapezoidal  integration  of  Equation  3  yielded  -S,  «  0.58  («  =  7,  95  percent  confidence  interval  0.53 

to  0.64)  and  S2  ~  1 .0.  The  shape  factors  are  hereafter  assumed  equal  to  unity  and  omitted  from  the 
equations  that  follow. 

The  conservation  of  sediment  flux  in  the  flow  is 


&.-BS 

dx 


(5) 


where  S  —  the  depositional  flux.  To  adjust  the  discharge  to  include  entrainment  and  deposition, 
such  that  deposition  does  not  decrease  concentration,  the  conservation  of  underflow  volume 
becomes 


^  =  B(E^.U-SIC,-PW,)  (6) 

where  =  the  entrainment  coefficient,  Q  is  the  deposit  solids  content,  and  P  is  the  depositional 
probability.  The  first  term  in  parentheses  on  the  right  side  of  Equation  6  increases  underflow  volume 
due  to  entrainment.  The  second  term  decreases  flow  due  to  deposition.  This  term  includes  both 
solid  and  liquid  phases  of  the  deposit.  The  third  term  in  parentheses  decreases  volume  due  to  liquid 
leaving  through  the  top  of  the  underflow  and  is  several  times  larger  than  the  deposit  volume  rate  of 
change.  The  underflow  is  not  an  irreversibly-mixed  solution  but  rather  a  settling  suspension.  By 
the  inclusion  of  the  second  and  third  right-hand  term.  Equation  6  is  different  from  equations  used 
in  other  models  of  particle-laden  gravity  currents  and  was  developed  under  the  assumption  that  the 
fluid  mud  layer  collapses  onto  the  bed  during  deposition.  Thus,  deposition  decreases  Q  (and  CQ) 
but  does  not  decrease  underflow  C  in  the  model. 

Wolanski  etal.  (1992)  showed  that  very  low  levels  of  turbulence  decrease  in  the  hindered  settling 
concentration  range  by  factors  of  2  to  1 0.  This  is  the  justification  for  the  inclusion  of  the  depositional 
probability  P  in  Equation  6. 

Equation  4  can  be  manipulated  to  give  the  thickness  change  along  jc  caused  by  momentum: 

—  -h  ®  I  (  (  Pi  Ap^<^(C<2) 

dx  [  h  h  [2  p  JQdx  [T  V  J  CQdx 

[2  JBdx]/^  ^ 
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where  Cy=  a  drag  coefficient,  the  factor  1.43  relates  the  bed  shear  stress  to  the  total  of  the  bed  and 
the  top  shear  stresses  (Findikakis  and  Law  1998),  and  the  bulk  Richardson  number  for  the  flow  is 

gAp/icos0  (g) 

pl/2 


The  excess  density  Ap  is  related  to  the  underflow  concentration  C;Ap  =  5'C  where 
=  (p^  -  pj  Yp^  and  p^  =  the  sediment  particle  density. 

Deposition.  The  depositional  flux  isS  =  where  the  probability  of  deposition  P  =  1  - 

Ted  =  the  threshold  shear  stress  for  deposition,  and  Q  =  the  concentration  at  the  interface  with  the 
deposited  bed.  Deposition  does  not  occur  at  bed  shear  stresses  greater  than  T^d  •  Since  the  imderflow 
was  found  to  be  highly  stratified  with  respect  to  C,  an  expression  from  Teeter  (1986)  was  used  to 
estimate  Cf 


Cu=C 


1  +  - 


Pe 


1.25  + 4.75  P 


,2.5 


where  Pe  =  W^h/K^  is  a  particle  Peclet  number,  and  = 0.067(Xf/p)^^2/i  is  taken  as  the  layer-average 
diffusivity.  The  valid  range  for  Equation  9  is  limited,  and  the  value  of  is  assumed  to  have  a 
maximum  of  4.  Since  underflow  concentrations  are  in  the  hindered  range,  Equation  7  decreased 
the  depositional  flux  relative  to  that  calculated  with  use  of  the  mean  concentration. 

Flow  regime.  A  Reynolds  number  (P)  criteria  for  the  turbulent-laminar  transition  has  been  proposed 
for  Bingham  plastic  materials  (Liu  and  Mei  1990),  and  found  to  be  applicable  to  mud  flows.  R  is 
composed  of  viscous  (Pp  and  yield-stress  (P.^)  components  depending  on  underflow  conditions, 
and 


P  = 


(l/P^,+l/P,) 


where  Pm  =  and  P^  = 


XyB^h^ 


(10) 


where  p  =  the  apparent  viscosity,  and  Ty  =  the  yield  stress.  Experimental  evidence  indicates  that  the 
turbulent-laminar  transition  occurs  at  P’s  of  about  2,000  (Liu  and  Mei  1990;  Van  Kessel  and 
Kranenburg  1996). 

Power  law  relationships  with  C  were  used  to  specify  p  and  x^ .  Van  Kessel  and  Kranenburg  (1996) 
and  Teeter  (1994  and  2000)  present  some  data  for  laboratory  clays  such  as  kaolinite  and  natural 
muds.  However,  data  for  the  range  of  C  measured  in  this  study  are  scarce.  Some  data  for 
low-concentration  fluid  muds  from  Gulfport  Harbor,  MS,  are  presented  in  Figure  15.  Those  data 
were  developed  from  shear-stress  sweeps,  starting  below  the  yield  stress,  with  a  Carri-med 
controlled-stress  rheometer  and  specially  collected  samples  (Teeter  1993).  A  concentric  cylinder 


19 


Figure  15.  Viscosity  at  50  sec'^  shear  rate  and  yield  stress  for  fluid  mud  from  Gulfport  Harbor,  MS 


geometry  was  used,  and  data  fit  to  a  Herschel-Bulkley  shear  stress  model  (Coussot  1994)  and  a 
Sisko  viscosity  model  (Barnes,  Hutton,  and  Walters  1989). 

The  yield  stresses  and  viscosities  at  50-sec''  shear  rate,  plotted  in  Figure  15,  are  higher  for  the 
overlapping  range  of  C  than  values  found  for  commonly  used  laboratory  clays  and  other  natural 
muds.  Measurements  of  yield  stress  were  obtained  for  three  channel  sediment  samples  from  the 
GIWW  about  6  km  north  of  Port  Isabel.  The  three  samples  were  from  the  top  7.5  cm  of  sediment 
cores  and  all  had  densities  of  about  400  kg/m^.  Yield  stresses  for  these  samples,  and  for  other  natural 
muds  including  Gulfport,  are  shown  in  Figure  16.  These  data  indicate  that  yield  stress  properties 
of  Laguna  Madre  mud  are  not  as  great  as  for  some  other  natural  muds. 

For  turbulent  conditions,  Cy  =  where  n  is  Manning’s  friction  factor.  For  laminar  conditions, 

the  drag  coefficient  was  estimated  as  the  maximum  of  the  turbulent  case  or  the  expression  presented 
by  Van  Kessel  and  Kranenburg  (1996)  for  laminar  flows: 


(11) 


Entrainment.  In  addition  to  deposition,  entrainment  of  overlying  water  can  reduce  C.  A  theoretical 
and  laboratory  investigation  of  entrainment  was  presented  by  Kranenburg  and  Winterwerp  (1997) 
and  the  resulting  relationship  between  entrainment  coefficient  and  Ri  tuned  for  fluid  mud  by 
Van  Kessel  and  Kranenburg  (1996)  is 
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Figure  16.  Yield  stress  for  Laguna  Madre  GIWW  sediments  and  other  muds 


where  Ri  is  defined  in.  Equation  9  and  the  condition  implies  that  turbulent  flow  is  required  for 
entrainment  to  proceed. 

Lateral  spreading.  A  buoyant-surface  spreading  rate  developed  by  Shirazi  and  Davis  (1974)  was 
adopted  to  the  case  of  negatively-buoyant  spreading  along  the  bottom,  yielding  the  expression: 


^  =  1.4 
dx 


B 


c^h  Ri 


-1/2 


(13) 


where  Cj  =  a  coefficient  inversely  related  to  the  concentration  difference.  An  upper  limit  of  2.0  was 
set  on  dB/dx  so  that  lateral  spreading  rate  could  not  be  larger  than  the  downslope  advance  rate. 

Initial  conditions.  A  flow  transition  formed  the  initial  condition  for  the  underflow  observed 
on  16  February.  A  horizontal  discharge  pipe  had  a  diameter  equal  to  an  appreciable  fraction  of  the 
receiving  water  depth.  A  turbulent  surface  flow  was  created,  which,  after  some  initial  entrainment 
and  spreading,  plunged  to  form  an  underflow  in  slightly  deeper  water  as  shown  schematically  in 
Figure  17  and  previously  in  the  Figure  8  photograph.  The  critical  Richardson  number  Ri  for  a 
plunging  underflow  described  earlier  is  slightly  greater  than  1  (Fang  and  Stefan  2000). 
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In  addition  to  the  pipeline  discharge  and  discharge  concentration  (Q,  and  C, ),  the  model  requires 
specification  of  the  initial  dilution  Sa,  initial  breadth  B^,  and  initial  Ri  which  were  estimated  from 
field  information.  The  initial  height  of  the  underflow  was  then  estimated  by 


Qo 


\2/3 


Pof?'o 
\S^0  j 


V/3 


(14) 


where  the  subscript  o  refers  to  values  at  the  underflow  transition,  =  SaQi,  and  Q  =  C,  /5a. 


Plunge  Point 


Figure  17.  Schematic  of  shallow  horizontal  pipeline  discharge,  turbulent  surface  flow,  and  transition  to 
underflow  at  plunge  point 


Model  Results.  Governing  and  auxiliary  equations  were  used  to  integrate  underflow  variables, 
starting  with  the  initial  conditions  described  in  the  last  section;  a  fourth-order  numerical  scheme 
was  used  to  calculate  Q,  CQ  ,  B,  and  h  at  1-m  intervals  along  the  downslope  trajectory  of  the 
underflow.  Depths  from  the  1-m  contour  toward  the  east-southeast  were  estimated  from  sample 
station  and  chart  depths  and  interpolated  out  to  a  distance  of 475  m.  Coefficients  c,  and  Manning’s 
n  were  taken  as  0.033  and  0.02,  and  model  results  were  insensitive  to  the  exact  values  used.  Example 
computed  profiles  of  underflow  and  deposit  heights  are  shown  in  Figure  18  for  1  and  6  hr  after  the 
discharge  began.  These  results  are  for  B^  =  45  m,  Sa  =  5.5,  initial  Ri  =  1.5,  =  0.5  m^/sec,  C)  = 

100  dry-kg/m3,  =  0.05  Pa,  and  the  settling  properties  described  previously.  The  rheological 
properties  were  set  between  those  of  kaolinite  and  Gulfport  Harbor  fluid  muds.  The  downslope 
extent  increased  between  hours  1  and  6  as  the  deposit  developed  and  reduced  slope  between  50- 
and  150-mx  distances.  Between  hours  1  and  6,  breadth- weighted  average  underflow  thickness  over 
the  footprint  decreased  from  0.78  to  0.61  m,  the  deposit  increased  from  0.05  to  0.28  m,  and  the 
footprint  area  increased  from  19,800  to  24,300  m^. 

The  underflow  extent  in  plan  view  is  shown  in  Figure  19  for  hour  6.  The  most  rapid  spreading 
occurred  as  Ri  increased.  Spreading  abruptly  stops  where  the  deposit  greatly  increases  bed  slope. 
The  variation  of  time,  velocity,  and  sediment  discharge  along  x  are  shown  in  Figure  20.  The 
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Figure  18.  Computed  underflow  and  deposit  profiles  at  two  times 

underflow  velocity  decreased  rapidly  near  the  origin,  then  more  slowly  over  the  remainder  of  the  x 
extent.  Deposition  starts  when  C/<  0. 1 0  m/sec.  Underflow  concentrations  were  1 8  dry-kg/m^  at  the 
origin  and  decreased  slightly  to  about  17  dry-kg/m^  at  the  maximum  downslope  extent  of  the 
underflow. 

ENTRAINMENT  INTO  WATER  COLUMN:  Entrainment  of  material  from  the  underflow  into 
the  water  column,  the  reverse  of  that  entrainment  previously  described,  was  implied  by  observed 
plumes  of  suspended  sediments  associated  with  the  location  of  the  underflow  footprint.  In  a  deeper 
estuarine  field  situation,  no  such  plume  was  observed  by  Thevenot,  Prickett,  and  Kraus  (1992). 
Though  currents  were  weak  at  the  Laguna  Madre  site,  winds  were  relatively  strong.  The  entrainment 
process  depends  on  the  local  momentum  balance,  turbulence  at  the  underflow  interface,  and  the 
magnitude  of  density  differences.  Data  collected  on  16  February  at  Laguna  Madre  were  used  to 
make  an  evaluation  of  the  entrainment  process  in  this  section. 

Details  of  the  underflow  interface  with  the  water  column  are  important  to  this  entrainment  process. 
Either  the  density  step  between  the  water  column  and  the  underflow,  or,  if  no  such  step  exists,  the 
gradient  at  the  top  of  the  underflow  is  used  to  scale  a  Richardson  number  (Turner  1986).  As 
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Figure  19.  Plan  view  of  computed  underflow  footprint  showing  spreading  with  distance  along  the  x-axis 


previously  discussed,  the  underflow  was  observed  to  be  highly  stratified  but  with  a  sharp  density 
jump  between  its  upper  surface  and  the  overlying  water  column.  Therefore,  the  magnitude  of  the 
jump,  along  with  the  depth  of  the  water  column  above  the  underflow  (/4),  was  used  to  scale  a 
Richardson  number.  The  magnitude  of  that  jump  was  much  smaller  than  the  overall  mean  density 
difference  between  the  two  layers. 

The  major  momentum  input  was  from  the  wind.  Wind  data  from  a  station  at  Rincon  maintained  by 
Texas  A&M  University,  Conrad  Blucher  Institute,  indicated  that  mean  wind  speeds  at  10-m  height 
(t4)  were  12.0  m/sec  for  1000  to  1200  CST  on  10  February,  and  were  7.6  m/sec  for  0900  to  1300 
CST  on  16  February  (standard  deviations  for  both  time  periods  were  about  1  m/sec).  An  approxi¬ 
mate  location  for  Rincon  is  given  in  Figure  1 . 

Hydrodynamic  forcing  was  assumed  to  equal  the  wind  stress,  some  of  which  goes  into  waves  and 
some  of  which  goes  into  currents.  Wind-waves  in  Laguna  Madre  tend  to  be  at  a  fully-developed, 
depth-limited  state  such  that  dissipation  was  nearly  equal  to  the  momentum  input.  Aalderink  et  al. 
(1985)  compared  two  models  which  used  wind  stress  directly  with  two  models  which  used  near-bed 
wave  orbital  shear  stress,  and  found  that  the  models  which  used  wind  stress  directly  were  better  able 
to  match  observed  TSM.  The  in-water  friction  velocity  was  estimated  from 

^*  =  {9aCdl9oi^lJa  (15) 


24 


ERDC  TN-D0ER-N16 
August  2002 


Figure  20.  Computed  variations  of  time  to  reach  x  distance,  underflow  velocity,  and  sediment  discharge 
along  the  x-axis 
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where  and  are  the  atmospheric  and  water  column  densities  and  Q  is  the  shallow-water 
atmospheric  drag  coefficient  taken  as  a  function  of  depth  and  wind  speed  (Teeter  et  al.  2001 ).  The 
average  value  for  on  16  February  was  about  0.0095  m/sec. 

At  high  interfacial  Richardson  numbers  (/?/*),  dimensionless  entrainment  {£)  is  the  result  of 
perturbations  in  the  interface  between  the  turbulent  water  column  and  the  underflow.  (Also 
assuming  that  the  molecular  Peclet  number  =  m,/,/v  is  greater  than  200  where  w,  and  /,  are  the 
turbulent  velocity  and  length  scales,  and  v  is  molecular  diffusivity.)  Under  conditions  of  turbulence 
vvithout  mean-flow,  the  laboratory  experiments  of  Long  (1975),  and  E  and  Hopfmger  (1986) 
confirmed  the  -3/2  power  law  described  by  Linden  (1973)  that 

E='^  =  KRifl^  (16) 

w* 


where  Ug  is  the  entrainment  velocity  or  the  downward  velocity  of  the  interface,  is  a  constant,  and 

the  interfacial  Richardson  number  is  defined  slightly  differently  from  Ri  as 


—  2 
P  M* 


(17) 


where  the  density  step  across  the  interface  Ap  =  (p  -  p  )/p ,  p  is  the  average  density  of  the  layers, 

and  Hg  is  the  depth  of  the  water  column  above  the  underflow.  I'he  scales  for  Ap  and  length  can  be 
chosen  differently  in  different  entrainment  systems.  Here,  although  the  underflow  is  stratified,  the 
mechanism  causing  that  stratification  involves  settling  and  not  diffusion  across  an  interface. 
Thicknesses  of  density  interfaces  are  typically  about  6  percent  of  the  depth  of  mixed  layers,  much 
thinner  than  the  stratified  underflow  layers  observed  here.  Values  of  Ri,  are  large,  and  interfacial 
perturbations  are  probably  intermittent,  consisting  of  vortex  rebounding.  Thus,  Ap  and  were 
scaled  by  the  overall  density  step  and  the  depth  of  the  water  column. 

Entrainment  and  deposition  to  the  underflow  by  settling  are  assumed  to  be  simultaneous  processes 
in  this  case.  Teeter  ( 1 994)  reviewed  laboratory  entrainment  experiments  involving  suspensions  and 
found  them  to  be  consistent  with  an  assumption  of  simultaneous  entrainment  and  settling.  Thus,  at 
a  depth-averaged  water-column  point  over  an  underflow 


H 


O 


-Fc 


(18) 


where  Q  is  the  depth-averaged  suspension  concentration  in  the  water  column,  t  is  time,  and  Fg  and 
Fg  are  the  entrainment  and  settling  flux  rates  at  the  interface. 

A  further  simplification  can  be  made  by  assuming  that  is  constant.  This  assumption  is  justified 
since  the  water-column  depth  is  small,  and  the  time  for  settling  or  turbulent  mixing  is  short  compared 
to  the  time-scales  for  underflow  spread  and/or  wind  speed  changes.  Under  equilibrium  conditions 
of  settling  and  entrainment,  Fg  and  F^  have  the  same  magnitude.  Furthermore,  the  settling  flux  for 
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the  water  suspension  can  be  estimated  by  use  of  laboratory  settling  tests  on  Laguna  Madre  channel 
sediments  described  by  Teeter  et  al.  (in  preparation)  and  shown  in  Figure  3.  These  results  were 
obtained  by  mixing  suspensions  with  site  water  and  allowing  them  to  settle  in  a  10-cm  diam  by 
1 .9-m-tall  column  under  quiescent  conditions.  Nine  initial  concentrations  were  tested,  and  settling 
velocity  (fV^)  was  found  to  increase  linearly  with  initial  concentration.  A  function  describing 
concentration-dependent  settling  rate  for  Laguna  Madre  channel  sediments  is 

=  ai  C",  0.1  <  Co  <  1  kg/ m^  (19) 


where  is  in  m/sec,  aj  =  0.806  x  10'^  m/sec,  and  the  exponent  n  =  L 

Settling  flux  depends  on  depositional  probability,  as  did  the  depositional  flux  S  described  in  the  last 
section.  Here,  turbulence  at  the  interface  is  assumed  to  be  low  and  intermittent,  and  the  depositional 
probability  is  assumed  to  be  unity.  The  depositional  flux  depends  on  the  near  bed  concentration 
Q,  but  in  this  case  it  is  assumed  that  Cf/C^  ~  1  •  Therefore,  Fg  =  UgC  and  Fg  =  W^Cg,;  a  simple 
model  for  the  water  column  suspension  concentration  is 


(k 

—  CUiiRi 


-3/2 

* 


l/(n+l) 


(20) 


Equation  20  was  recast  to  solve  for  the  entrainment  coefficient  K  with  field  data.  The  Laguna  Madre 
measurements  and  Equation  7  indicated  that  C/Q  =  1.07,  and  column  Q  values  were  adjusted 
accordingly.  In- water  friction- velocities  were  calculated  using  Equation  19  and  a  constant  wind 
speed  of  7.6  m/sec.  Underflow  concentration  at  the  interface  was  assumed  to  be  3  dry-kg/m^. 
Results  for  K  and  other  select  parameters  are  presented  in  Table  2. 

The  flux  Richardson  number  (RiJ)  was  calculated  as  {uJu^)Ri^  (Turner  1986)  and  represents  the 
ratio  of  turbulent  kinetic  energy  dissipation  by  buoyancy  flux  and  turbulent  production.  Flux 
Richardson  numbers  greater  than  0. 1  are  associated  with  damping  of  turbulence,  and  the  magnitudes 
of  iJiyobtained  here  indicate  appreciable  suppression  of  turbulence  at  the  interface. 

The  K  (and  TSM)  values  were  log-normally  distributed.  The  median  K  value  was  2.8  in  fair 
agreement  with  the  laboratory  result  of  3.8  reported  by  E  and  Hopfinger  (1986).  Based  on  this 
estimate  for  K,  an  estimate  of  Q  was  made  for  1 0  February  with  Equation  20.  Assuming  underflow 
conditions  were  the  same  as  before,  where  Ap/p  =  0.00081, ,  and  that  Hg  =  2m  and  =  0.0174 
m/sec,  then  the  Ri^,  for  that  day  was  about  52.  Equation  20  predicts  Q  =  690  mg/L  or  somewhat 
higher  than  the  high  values  observed  in  the  field  (Figure  5). 

DISCUSSION  AND  CONCLUSIONS:  Some  of  the  important  results  from  the  field  measure¬ 
ments  were  as  follows:  the  underflow  thickness  was  relatively  uniform  and  decreased  rapidly  near 
the  limit  of  downslope  extent.  The  underflow  had  a  distinct  upper  surface  or  interface  with  the 
ambient  water  column,  but  this  concentration  was  only  about  3  dry-kg/m^.  Underflow  layers  were 
sediment  and  density  stratified.  Turbulence  in  the  underflow  was  not  sufficient  to  mix  sediment 
vertically.  Thick  deposits  formed  under  the  underflow.  The  implication  of  these  observations  was 
that  the  underflow  was  slow  moving.  An  apparent  absence  of  appreciable  entrainment  also  indicated 
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Table  2 

Entrainment  Conditions  at  Underflow  Profile  Stations 


Time,  CST 

. 

TSM,  mg/L 

Ho,  m 

Ri. 

Fsx  10= 

kg/m^/sec 

K 

Rif 

1057 

168 

0.86 

86 

2.6 

0.7 

-0.,  1 

1110 

240 

0.26 

133 

5.3 

2.8 

0.2  1 

1122 

308 

0.41 

69 

8.8 

1.8 

0.2 

1136  I 

228 

0.37 

128 

4.8 

2.4 

0.2 

1216 

- ] 

1  480 

0.29 

134 

21.3 

11.4 

1 

1233 

708 

^  0.2 

124 

46.3 

22.1 

2  I 

1308 

260 

0.25 

151 

6.2 

4 

0.3  J 

slow  underflow  movement  and  high  Richardson  numbers.  The  concentrations  near  this  upper 
surface  and  layer  average  values  were  relatively  uniform  along  the  length  of  the  underflow.  This 
observation  indicated  the  underflow  collapsed  vertically  while  depositing  and  led  to  the  model 
feature  that  tended  to  maintain  underflow  concentrations  while  deposition  occurred. 

The  results  from  the  model  were  in  general  agreement  with  features  observed  in  the  field.  The  model 
predicted  the  underflow  to  be  slow  moving,  starting  turbulent  and  becoming  laminar  after  85  to  150  m 
distance  downslope.  Calculated  underflow  concentrations  were  almost  constant.  Deposition  started 
near  the  source,  became  appreciable  after  100  m,  and  decreased  rapidly  at  the  downslope  extent  of 
the  underflow.  Underflow  thicknesses  in  the  model  were  similar  to  those  observed  in  the  field.  The 
extent  of  the  underflow  could  not  be  accurately  measured  in  the  field  but  appeared  to  be  in  rough 
agreement  with  the  model  predictions.  Lateral  spreading,  which  was  controlled  by  bed  topography, 
was  captured  only  in  a  very  approximate  way  in  the  model.  The  model  did  not  include  path 
switching  that  apparently  occurred  in  the  field  nor  the  complexin^  of  bottom  topography  that 
influences  underflow  spreading. 

Laboratory-determined  settling  velocities  were  used  in  the  description  of  imderflow  spreading  and 
in  the  analysis  of  water  column  entrainment  of  underflow  material.  Laboratory  measured  values 
are  probably  not  the  same  as  field  values,  but  obtaining  measurements  in  the  field  is  problematic. 
Teeter  (2001)  showed  that  for  low-concentration  Laguna  Madre  suspended  sediments,  quiescent 
column  tests  yielded  settling  rates  representative  of  disrupted  floes.  Very  mild  turbulence  produced 
much  larger  values,  but  shear  rates  greater  than  2  sec'*  produced  settling  rates  that  were  not  much 
different  from  quiescent  values.  The  uncertainty  in  fV^  affects  the  results  of  both  the  underflow  and 
water  column  analyses,  if  performed  in  a  predictive  mode. 

Our  understanding  of  fluid-mud  flow  properties  is  incomplete  and  measurements  are  difficult.  The 
existence  of  a  yield  stress  may  lead  to  an  unsheared  plug  flow  zone  (Coussot  1 994)  in  the  underflow 
and  could  make  R  values  based  on  layer  average  properties  unrepresentative  of  interface  conditions 
where  entrainment  occurs.  In  the  present  case,  the  underflow  layer  was  highly  stratified,  and  the 
velocity  profile  was  difficult  to  evaluate.  The  flow-regime  evaluation  imposed  on  Equation  1 2  only 
captures  the  effect  of  fluid-mud  flow  properties  on  entrainment,  if  such  properties  are  known  or  can 
be  estimated. 
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Rheological  data  on  muds  are  relatively  scarce.  While  varying  rheological  parameters  over  an  order 
of  magnitude  from  those  used  in  the  example  calculation  did  not  greatly  affect  results,  increases  in 
viscous  and  yield  stress  properties  eventually  caused  the  model  underflow  to  freeze  and  become 
numerically  unstable.  The  range  from  laboratory  clays  to  cohesive  natural  muds  such  as  those  in 
Figures  15  and  16  affect  computed  underflow  R,  friction,  flow,  and  ultimate  extent  of  spreading. 
The  predictive  capability  of  the  model  is  therefore  dependent  on  rather  extensive  site-specific  field 
information.  Important  factors  include  sediment  composition,  settling  and  rheological  charac¬ 
teristics,  bed  topography,  ambient  currents,  winds,  and  waves. 

Elevated  water-column  suspended  sediment  concentrations  were  caused  by  underflow  entrainment 
into  the  water  column  by  wind-wave  forcing.  Entrainment  model  coefficients  were  consistent  with 
previously  reported  values  for  high  Ri^  situations  when  used  with  wind-stress  forcing. 

The  pipeline  discharge  underflow  represents  the  greatest  potential  for  local  turbidity  generation,  if 
it  is  entrained  into  the  overlying  flow,  since  it  contains  the  vast  majority  of  sediment  particles 
discharged.  Field  observations  indicate  that  at  times  of  high  bed  shear-stress,  entrainment  of 
underflow  material  can  generate  a  turbid  plume  extending  some  distance  from  the  discharge,  but 
not  necessarily  downstream  from  the  discharge.  Thus,  the  area  of  concern  with  respect  to  water 
column  impacts  of  a  pipeline  discharge  is  not  confined  to  the  vicinity  of  discharge,  but  also  includes 
the  area  over  the  underflow  that  might  extend  hundreds  of  meters. 
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